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ABSTRACT 

Context. Fomalhaut is a young (2 ± 1 x 10** years), nearby (7.7 pc), 2 M© star that is suspected to harbor an infant planetary system, interspersed 
with one or more belts of dusty debris. 

CNj Aims. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7" and 36.7" at wave- 
lengths between 70 fim and 500 yum. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. 

' . ' , The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the 
star, excess emission is detected. We aim to construct a consistent image of the Fomalhaut system. 

Methods. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models 
of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to 
the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be 
blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. 
^ Results. The appearance of the belt points towards a remarkably active system in which dust grains are produced at a very high rate by a collisional 
$^ _ cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. 
, The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 1 10 Earth masses. From compar- 
' ison of their scattering and thermal properties, we find evidence that the dust grains are flufi"y aggregates, which indicates a cometary origin. The 
I I excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free 
emission from a stellar wind. 

^ Key words. Stars: individual: Fomalhaut; Circumstellar matter; Planetary systems; Radiative transfer; Zodiacal dust 
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O 1. Introduction uncertain jKalas et al.ll2008ll201 UUanson et al.ll2012h . The most 
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prominent feature in the disk is the main belt, which is located 



^ Fomalhaut was one of the first mam-sequence stars shown to o^tsj^e of 130 AU. The center of this belt is shifted with respect 

O have aninfrared excess caused by dust grams m orbit around the j^e star, implying significant eccentricity, probabl y caused by 

(N| . star (| Aumann| |1985). Further studies of this phenomenon have gravitationa l interaction w ith one or more planets kalas et alJ 

T— I I shown that approximately 15% of s olar-hke main-seque nce stars |2005ll2008tlChiang et al Il2009h 
L| . have a detectable amount of debris (iBrvden et al.l2006l) . Belts of 
. !^ ' debris typically reside in planetary systems at locations where no 

planets have formed because either the formation time scale for 2. Observations 

$-H planets was too long, such as in the Solar System's Kuiper Belt, o i u i i ■ 

. or because other planets in the system have stirred up the belt Herschel images 

before a plane t could be formed, such a s in the Solar System's Far-infr ared images of Foma lhaut have bee n obtained with the 

Asteroid Belt (|Wyatt||2008HKrivov||2010')- PACS (PogHtsch etaP 12010 ) and SPIRE dGriffin et all 120101: 

Fomalhaut's disk was previously imaged from the ground at Swinvard et al ,20Tol) instruments aboard the Herschel Space 



450 fim and 850 fim, and from space with the Spitzer and the Observatory (iPilbratt et al.l 1201 Ol) in five photometric filters 



Hubble Space Telescopes (iHoUand et al.ll2003HStapelfeldt et al] centered around 70 yum, 160 yum (PACS), 250 yum, 350 jim 

[2004; Kalas et al. 200^. A planet candidate has been reported and 500 //m (SPIRE). The Observation IDs are 1342211953, 

in the disk, though 3'^'* epoch confirmation of this object is still 134221 1954 (PACS) and 1342195939 (SPIRE). The full width 

at half maximum of the point spread functions (PSFs) indicate 

* Herschel is an ESA space observatory with science instruments a spatial resolution of the obsei-vations of 5.7", 11.3", 18.1", 

provided by European-led Principal Investigator consortia and with im- 25. 1" and 36.7" for the five bands respectively. This coiTesponds 

portant participation from NASA. to 43, 86, 140, 190 and 280 AU at the distance to Fomalhaut. 
** Postdoctoral Fellow of the Fund for Scientific Research, Flanders The PACS data were proces sed in the Herschel Interactive 

*** Aspirant Fellow of the Fund for Scientific Research, Flanders Processing Environment (HIPE, IOttll201Cil) version 6.0, apply- 
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Fig. 1. Herschel PACS 70 /vm image of Fomalhaut. The positions 
of the star and of the planet candidate Fomalhaut b are indicated. 
The hatched circle indicates the spatial resolution (5.7"). The 
black-blue-white color scale is linear. North is up, east to the 
left. 



ing the standard pipeline steps. The flux conversion was done 
using version 5 of the response calibration. Signal glitches due 
to cosmic-ray impacts were masked using the PACS photMMT- 
Deglitching task in HIPE on the detector timeline. Then a first 
coarse map was projected. A high-pass filter was applied to re- 
move the low frequency drifts. To avoid the introduction of arte- 
facts, we used the coarse map to mask out the disk from the de- 
tector timeline prior to the application of the filter The scale of 
the high-pass filter was 15 frames for the blue maps (70 //m), and 
25 frames for the red maps (160 yum). The detector time-series 
signals were then summed up into a map using the PACS phot- 
Project task. The pixel scale for the 70 /vm map was set to 1", 
while the scale for the 160 jim map was 2". The source was cov- 
ered in two independent observations, scanning the sky in two 
orthogonal directions. We combined the two detector time series 
and projected these together into the final maps. 

Also the SPIRE data were reduced using HIPE. The high 
signal-to-noise ratio of these data allowed us to created over- 
sampled maps with pixel scales of 3", 4.2" and 7" in the 250 fim, 
350 /vm and 500 /vm bands respectively, equating to approx- 
imately 6 samples per beam FWHM. The standard pipeline 
script provided with HIPE was used to reduce the data, with the 
naiveMapper task being used to construct the maps. A simple 
linear baseline subtraction was the only filtering performed on 
these data. 

Table [1] summarizes the pixel scale, total map integration 
time and signal-to-noise ratio of the projected Herschel images. 
We also indicate the integration time per pixel, which is the total 
time a pixel in the central part of the projected map was seen by 
a detector pixel. 

Figs. [T] |2] and [3] show that the system is well resolved at 
70 ;um, barely resolved at 500 jim, and marginally resolved at 
intermediate wavelengths. A first qualitative look at the images 



shows that the thermal emission of the dust in the belt around 
Fomalhaut appears smooth and not clumpy, contrary to what 
was expected by Wvatt & Dent (2002). The smoothness hints at 
a high dust replenishment rate by numerous collisions, rather 
than sporadic replenishment by rare collisional events. The ap- 
parent azimuthal brightness variations are due to geometric pro- 
jection effects. Moreover, distances from the belt to the star vary 
along the belt, as it is eccentric. The southern ansa is located 
closer to the star than the northern ansa, and is therefore warmer 
and brighter This is the so-called pericente r glow, first ob- 
served in a debris disk bv lTelesco et al.l (l2000l) . The planet can- 
didate Fomalhaut b may produce local density enhancements. 
However, there is no evidence of an excess or a lack of dust 
emission near the position of Fomalhaut b, nor at any other lo- 
cation along the belt. 

In the 70 jdva image, the star shows up as an unresolved point 
source. The flux of this central source, however, exceeds the stel- 
lar photospheric flux. This is in agreement with excess flux close 
to the star previousl y found at near-inf rared wavelengths and at- 
tributed to hot dust (lAbsil et al.ll2009h . The region between the 
unresolved central source and the belt is not devoid of dust ei- 
ther; the 70 yum image shows emission in this inner disk as well. 
The dust particles in this region either originated in the belt and 
spiraled in past the orbit of Fomalhaut b, or were released from 
an unresolved inner belt and moved outwards. 

2.2. Belt geometry 

The geometric properties of the belt were derived from the PACS 
70 jum and 160 images in the following way. Along each 
pixel row and column intersecting the belt, the coordinates of 
the pixel with the highest surface brightness are determined. In 
this way, the position of the ring in the image is captured by a 
set of pixel coordinates {x,y) along the ring. An ellipse is fitted 
to these coordinates, which yields the major axis, the minor axis, 
and the position angle of the belt as it appears on the sky. The 
position of the star, relative to the belt center, is measured from 
the PACS 70 ;um image by fitting a two-dimensional Gaussian 
curve to the central unresolved point source. From the parame- 
ters of the apparent ellipse, the stellar position, and the known 
distance to Fomalhaut, the orbital elements of the belt are com- 
puted: the semi-major axis a, the eccentricity e, the inclination 
/, the longitude of the ascending node Q., and the argument of 
periastron o). 

We repeat this ellipse-fitting procedure 1000 times by adding 
random noise in each spatial pixel (spaxel) of the image. The 
magnitude of this added noise is equal to three times the pixel- 
to-pixel noise of the original image, estimated by taking the stan- 
dard deviation of the image outside the location of the Fomalhaut 
system on the sky. As a result of this added noise, the set of belt 
locations {x,y) are not exactly the same for two dififerent sim- 
ulations (although several pixels may be in common). Also the 
measured stellar position varies. Each simulation will therefore 
yield slightly different parameters. Our final set of belt parame- 
ters are the mean values, and the errors the standard deviations, 
of the outcome of these 1000 simulations. The same procedure 
is applied to derive the belt parameters from the PACS 160 jjm 
image. Because the star is not visible in this image, however, we 
use the relative stellar position measured in the 70 yum image. 

The geometric properties of the belt are summarized in 
Table |2] They are similar to those derived from optical scattering 
images (Kalas et al. 2005), apart from the argument of periastron 
(jj. We find that the star is located almost perfectly on the appar- 
ent major axis of the belt on the sky, leading to w = 1° + 6°. 
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Table 1. Herschel images. 



Image 


Pixel" (") 


Integration time (s) 


SB" 


SNR" 






Total map 


Per pixel*^ 


(mJy/"2) 




PACS 70 jum 


1 


10956 


15.6 


16.5 


77 


PACS leOyum 


2 


10956 


16.7 


6.28 


59 


SPIRE 250 yum 


3 


2906 


2.1 


1.90 


32 


SPIRE 350 yum 


4.2 


2906 


2.8 


0.59 


21 


SPIRE 500 yum 


7 


2906 


3.6 


0.19 


27 



" Pixel size of the projected image. 

" Surface brightness and signal-to-noise ratio in the brightest pixel of the image. 
° In the central part of the image where the source is located. 



Table 2. Geometric properties of the elliptical belt. 





Optical" 


70 yum 


160 //m 


Semi-major axis a (AU)'' 


141 +2 


137.5 ±0.9 


125 ± 3*^ 


Eccentricity e 


0.11 ±0.01 


0.125 ±0.006 


0.17 ±0.07 


Inclination / (°) 


65.6 ± 0.4 


65.6 ± 0.5 


64 ±2 


Longitude of ascending node (°) 


156.0 ±0.3 


156.9 ±0.5 


160 ±2 


Argument of periastron 


31 ±6 


1 ±6 


d 


Offset between star and ring center (AU)'' 


15 ± 1 


17.2 ±0.9 


d 



" iKalas et alj f2003) . 

^ At peak surface brightness. 

Spatial resolution of the 160 yum image is insufficient to derive a reliable semi-major axis. 
'' Star not visible at 160 /jm. 

In the plane of the belt. 



iKalas et al] j2005t) do not explicitly mention the value, but w 
can be estimated from the measured difference (6 - 14°) be- 
tween their longitude of the nodal line (Q - 156°) and their belt 
center-star offset direction (170°), deprojected for the inclina- 
tion (i = 65.6°): CO - arctan(tan5 /cos;) = 31 ± 6°. The reason 
for this discrepancy is unclear to us, but we are confident in our 
value. Our Herschel 70 jim image is unique in that it captures 
both the stellar light and the thermal dust emission simultane- 
ously, in contrast to the Hubble Space Telescope (HST) images, 
where the direct starlight had to be blocked to be able to detect 
the light scattered by the belt. Our measurement of w is there- 
fore more direct. To elaborate on this issue, we have tested both 
values of co with our radiative-transfer model (see Sect. |3). In 
the case of w - 31°, the pericenter glow in the model images is 
clearly Oj^the major axis of the belt, towards the South, unlike 
the observed on-axis pericenter glowQ. The on-axis position of 
both star and pericenter glow, lead us to conclude that co must be 
close to zero degrees. 

Although the HST optical images have a spatial resolution 
much higher than that of Herschel/PACS at 70 fim (0.5 AU vs 
43 AU at the distance of Fomalhaut), the error bars on the belt 
parameters are comparable. This is mainly due to the superb 
signal-to-noise ratio of the PACS images (see Table[TJ. 

The belt properties derived from the 160 /urn image are less 
certain. The derived semi-major axis (125 AU) is significantly 
smaller than the one derived from the 70 /im image (137.5 AU). 
We ascribe the difference to the lower angular resolution; The 
larger PSF at 160 /urn makes the belt appear smaller, because it 
squeezes the locations of peak surface brightness inwards. 



' The model images are convolved with the PACS 70 yum point spread 
function before comparison with the observations. All asymmetries re- 
lated to the latter are hence taken into account. 
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Fig. 2, Cut through the HerschelfPACS images along the belt's 
major axis on the sky. The surface brightness profiles are normal- 
ized by the maximum values indicated in Table[T] The off-center 
position of the star is indicated. The belt's peak positions are lo- 
cated closer to the belt center at 160 //m, which is a consequence 
of the reduced angular resolution (see Sect. 12. 2t . 

2.3. Spectral energy distribution 

Herschel photometry was derived from the images described 
above, by integrating the flux within a circular aperture around 
the geometric center of the belt on the sky. We used an aper- 
ture radius of 52", i.e. 400 AU at the distance of Fomalhaut 
(Table |3]l. The mentioned uncertainties consist of the statistical 
errors, amounting to ~10% of the measured flux. According to 
the SPIRE and PACS documentation available from the Herschel 
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Fig. 3. Herschel far-infrared images of Fomalhaut. (a) Observations, (b) Model images, (c) Absolute value of the residuals. The 
hatched circles indicate the spatial resolution. The color scale is linear rainbow to increase contrast at low surface brightness. The 
scale ranges from zero to maximum surface brightness (see TablelTJ of the image in columns (a) and (b); In column (c), the color 
scale is the same as in column (a). North is up, east to the left. 
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Table 3. Herschel photometry of Fomalhaut. 



Band 


Intej 


;rated 




Flux contribution (%)'' 






flux 


(Jy)" 


Point source 


Inner disk 


Belt Blow-out grains 


PACS 70 fim 


10.8 


+ 0.9'^ 


5 


21 


74 


27 


PACS 160 yum 


6.2 ± 0.6 


2 


23 


75 


7 


SPIRE 250 yum 


2.7 d 


:0.3 


2 


25 


73 


3 


SPIRE 350 yum 


1.1 d 


:0.1 


2 


27 


71 


2 


SPIRE 500 yum 


0.50 


±0.05 


2 


29 


69 


2 



" Within a circular aperture with a radius of 52" (400 AU) around the geometric center of the belt. 
From the best-fit model. 

The excess flux at the location of the star is 50 + 10% of the stellar flux at 70 yum, or 0. 17 ± 0.02 Jy. 



Science Center web pagefl the absolute calibration error on the 
photometry is another 10%. 

Additional photometry of the Fomalhaut system was 
extracted from several online catalogues, available from 
|http://vi z ier.ustrasbg.fr/viz-bin/Viz ieR Stellar photosphere 
models (ICastelU & Kuruczl l2003h were fitted to the optical 
Geneva photometry ( Rufenerl " IQSSV We have interpolated 
between the grid models of [Castelli & Kurucz (200 31), ac- 
cording to the method explained in IDegroote et al.l (1201 ll 
their Sect. 2). The derived fundamental parameters of the star 
are an effective temperature Teff = 8600 + 200 K, a surface 
gravity log^ = 4.1 + 0.2 dex, and an interstellar reddening 
E(B-V) < 0.08. The photometric angular diameter of the 
star is 2.20 + 0.04 milli-arcsec , in perfect ag r eemen t with 
interferometric measurements (iDi Folco et al.l l2004l) . The 
infrared photometry is taken from the A KARI/IRC and FIS 
All-Sky Survey point source catalogues (llshihara et al.l I20T0I: 
[Yamamura et al. 2010 ), the ESO catalogue IR photometry of 
calibrator stars (Ivan der Bliek et al. 1996), and the Spitzer Space 
Telescope MIPS catalogue (Su et al., .2006) . The two sub-mm 
measurements, at 450 yum and_850 um, are consistent with the 
SPIRE photometry OHoUand et al.ll2003h . The spectral energy 
distribution (SED) of the Fomalhaut system is shown in Fig.jS] 



3. Modeling 

3.1. Three components 

We use a three-component model to reproduce the detected ther- 
mal dust emission. A central, unresolved component responsible 
for the stellar plus excess emission at the location of the star, a 
narrow ring of dust produced in a collisional cascade for the belt, 
and a power-law surface density to reproduce residual emission 
inside the belt. 

We assume that the dust in the belt originates from a narrow 
source region, a ring between radii ri and ri in which dust grains 
are produced through collisions with a size distribution accord- 
ing to the equilibrium cascade: f{s) oc s^^-^, with s defined as the 
radius of the grain (Dohnanyi 1969). The particles are assumed 
to be released from a parent body on an orbit inside the ring. 
For small grains, their orbits are modified by stellar radiation 
forces on the particles. The modification depends on the ratio 
P - F^idlFgaw of radiation and gravitational forces acting on 
the particles. Particles with fi < 0.5 will remain in bound orbits 
around the star, with periastrons equal to their release points 
in the ring, and forced eccentricities e - /3/(l - /3) so that the 
semi-major axes become a - rp/(l - e). Hence, particles with 



http://herschel.esac.esa.int/ 



increasing /3 will find themselves in orbits with increasing ec- 
centricities and apastron distances. Particles with a given will 
produce a surface density distribution inversely proportional to 
2:7r^fl -sja^e^ - (r - op- in the region a{\ — e) < r < a{\ + e). 

Particles with fi > 0.5 will be on hyperbolic orbits, and grad- 
ually leave the system. These so-called blow-out particles are 
constantly replenished by the collisional cascade and contribute 
at low levels to the surface density outside of the ring. We as- 
sume that the small-grain extension of the collisional cascade 
power-law is replenished every fiepi years, which provides a scal- 
ing factor for the surface density of the blow-out particles. The 
replenishment time is a free parameter of the model. 

For particles in the inner-disk region, we take a power-law 
size distribution with f{s) oc s^^. Such a distribution will result 
from drawing grains by means of Poynting-Robertson (PR) drag 
from the belt (Burns et al. 1979), as explained in Appendix lAl 
Moreover, we remove all particles with p > 0.5 from the inner 
disk, since these will be blown out quickly. The size distribu- 
tion in the inner disk is not well constrained by the Herschel 
images, however. Besides PR drag, evaporating comets and ad- 
ditional debris belts closer to the star are possible sources of dust 
grains in the inner disk between star and outer belt that could 
have very different size distributions. Including only PR drag 
and assuming that no collisions occur in the PR stream inward 
of the source ring, a constant surface density in the inner disk 
is expected (see Appendix lAl). resulting in a surface brightness 
profile that is heavily peaked towards the central star but this 
is inconsistent with our observations. A planetary system inside 
the ring could clear out the particles as they drift in, creating a 
surface density profile that decreases towards the star. We model 
this distribution using a power-law surface density where the ex- 
ponent is a free parameter. 

The central source in the 70 yum image is modeled with a 
point source, convolved with the PSF of the PACS instrument 
at this wavelength. The model hence consists of three geometric 
components: a belt, an inner disk and a central point source. 

3.2. The dust model 

The dust particles in our simulations are mixtures of common 
astronomical species, including water ice, following the compo- 
sition proposed by Min et al. (20ll]). This is a mixture consisten t 
with the solar abundances derived bv Greve sse & Sauvail(ll998h . 
The resulting composition contains by mass 32% silicates, 10% 
iron sulfide, 13% amorphous carbon, and 45% water ice. 

To simulate inhomogeneous grains we mix the various ma- 
terials using effective medium theory. F or this we use the 
Bruggeman mixing rule (Bruggeman 1935: iVoshchinnikov et al.l 
120071) . We simulate porosity of the grains by adding a fraction of 
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Fig. 4. The mass absorption coefficients of the dust grain mixture 
used in the models, averaged for grains with sizes below (full 
line) and above (dashed line) the blow-out size (13 /vm). See text 
for details. 

vacuum to the mixture. The grai n shape is descr ibed by a distri- 
bution of hollow spheres (DHS. lMin et al.ll2005 ) with an irreg- 
ularity parameter /n^x = 0.8. With a mild porosity of 25% vac- 
uum, we find that the blow-out size for grains in the Fomalhaut 
system is 13yum. We tested various values of the porosity and 
also d ust particles without ice to mimic the effect of photosput- 
tering (i Grigorieva et al.ll2b07 !). 

For the initial size distribution formed by the collisional cas- 
cade we use a power law, as stated above, a minimum grain size 
of 1 /im, and a nominal maximum grain size of 5000 /vm. The 
latter value is not well constrained by the observations. The size 
distribution at each location in the disk is computed following 
the dynamical model described above. The averaged mass ab- 
sorption coefficients for grains formed in the collisional cascade 
with sizes below 13/im, i.e. the blow-out grains, and the aver- 
aged opacities for grains larger than the blow-out size are plotted 
in Fig.|4| Note that these opacities do not necessarily represent 
the averaged opacities at any location in the disk since the size 
distribution is significantly altered by the radiation dynamics. 

3.3. The model grid 

We set up a grid of models with varying inner and outer radius 
for the ring where the planetesimals collide. We fixed the value 
of ri equal to the inner-edge radius of 133 AU that was well de- 
termined based on the HST images (Kalas et al. 2005). We also 
tested a slightly different ri of 123 and 143 AU. For the outer 
radius of the ring, we took r2 = ri 4- [2,5, 10,20,30,40] AU. 
For the inner disk, we tested surface-density power laws propor- 
tional to r" (PR), r' and r^. The outer edge of the inner disk is 
set to the inner radius of the ring. On the inside, we assume that 
our icy dust particles sublimate when they reach the sublimation 
temperature of water ice, i.e. 100 K in vacuum. This effectively 
cuts off the inner disk at an inner radius of ~35 AU. 

3.4. Fitting approach 

Model images in each of the five Herschel bands were created 
using the radiative-transfer code MCMax (Min et al. 2009|). The 
code was modified to allow for an eccentric belt, with the star 
in one of the foci. The model images are convolved with the 



measured PSF of the PACS and SPIRE instruments at each ob- 
servational wavelength. PACS and SPIRE images are obtained 
in broad-band filters. It is therefore not correct to directly com- 
pare the monochromatic model images with the observations. 
We take the filter transmission into account by computing model 
images at three wavelengths within each of the PACS and SPIRE 
bandfl The three images, weighted with the filter transmission 
curve, are then combined into a single image. The combined im- 
age represents the band-integrated model image. 

The images of the three geometric components are stored 
separately. Hence, the images of the belt, the inner disk and the 
unresolved point source of a given model can be independently 
scaled to the observations. We start with the PACS 70 fim image 
and determine the best-fit linear combination of the three com- 
ponents using a non-negative least-squares (NNLS) fitting rou- 
tine (Lawson & Hanson 1974). Because of the detected excess 
flux close to the star, the flux of the central source is not fixed 
to the stellar flux, but is allowed to exceed it. The circumstel- 
lar material around Fomalhaut is optically thin at all considered 
wavelengths. The scaling of the inner disk and belt is therefore 
equivalent to scaling the total dust mass in these components. 
The linear combination of the three geometric components hence 
provides the best possible fit to the 70 fim image, for each model 
in the grid. 

Two quantities are computed to assess the goodness-of-fit. 
First, a Pearson chi-square, summed over all spatial pixels where 
the surface brightness is larger than 20cr, judges the general ac- 
cordance between the linearly scaled model and the observation. 
Second, the maximum of the norm of the residuals is computed, 
to control the most extreme deviation between model and obser- 
vation. Similar values were determined for the Herschel images 
at longer wavelengths. While these images are not used in the 
NNLS fitting routine, the agreement between model and obser- 
vations is checked a posteriori and taken into account to deter- 
mine the best-fit model in the grid. 

3.5. Best modei 

The best-fit model has a Pearson chi-square value, summed over 
all five images, of 0.42. The maximum pixel-to-pixel residual 
occurs in the 250 fim image, where it amounts to 24% of the 
maximum surface brightness of the observation. The model im- 
ages and residuals are shown in Fig. |3] The general accordance 
between model and observations is very good. The most signif- 
icant residual structure is seen in the 250 fim image, where the 
belt appears to be more extended than the model predicts. 

The best model has a debris source ring located between ri 
= 133 AU and r2 - 153 AU. The belt contains a total mass of 8 x 
10^^ g in dust grains up to sizes of 5000 fim. The replenishment 
time, determining the amount of blow-out grains present in the 
system, is very short, approximately 1700 yr, leading to a total 
mass of 3 X 10^"* g in grains with sizes below 13 fim. The inner 
disk contains 4 x 10^^ g of dust grains, and is responsible for 
21-29% of the flux in the Herschel images (TableO. Its surface 
density increases outward, linear with r, to the inner radius of the 
ring. The central unresolved excess-emission component has an 
integrated flux that amounts to 50% of that of the star at 70 jum 
(TableO. 

The grain properties that best fit the SED and images are 
icy grains with a vacuum fraction of 25% by volume. However, 



^ Three wavelengths per band is sufficient to sample the smooth 
wavelength-dependence of the model within the band, while keeping 
the computational time limited. 
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Fig. 5. The infrared SED of the Fomalhaut system. Black dots 
are photometric measurements from the literature, red dots the 
Herschel photometry. The full red line represents the best model, 
the dotted and dashed red lines indicate the contribution of the 
large grains and blow-out grains respectively. The blue line is 
the stellar photosphere model. The green dot indicates the total 
flux of the central point source in the 70 /im image, and shows 
the excess of the stellar photosphere. The excess-to-stellar lumi- 
nosity ratio, estimated from the model, is 7.4 x 10"^, of which 
27% is produced by blow-out grains. 



we find Uttle difference in total dust mass in blow-out grains - 
and hence in replenishment time- between models with different 
grain properties. This was to be expected; The key physical pa- 
rameter that determines the thermal and dynamical properties 
of the dust particle is not the size, but the absorption efficiency. 
The grain model merely couples the absorption efficiency to a 
specific grain size, which is different for e.g. a different degree 
of porosity. Although grain models link different grain sizes to a 
similar absorption efficiency, it is this latter quantity that deter- 
mines the appearance of the model. We conclude there is a sig- 
nificant amount of mass contained in grains with sizes smaller 
than the blow-out size in the Fomalhaut system. This conclusion 
is robust against assumptions for the grain structure and compo- 
sition. 



4. Results 

4. 1 . Hot excess emission 

Excess emission clos e to the star has been reported at near- 
infrared wavelengths (l Absil et ani2009l) . The SED of our model 
shows that the emission from the inner disk and belt only be- 
comes visible at wavelengths beyond 20 fim. Closer inspection 
of the literature photometry shows that there is indeed a contin- 
uous excess above the stellar photospheric flux in the near- to 
mid-infrared wavelength range. In Fig. |6] the SED of this excess 
is shown. The PACS 70 jim photometric data point of the on-star 
excess is added. The error bars include the measurement uncer- 
tainties as listed in the different catalogues, but also the spread 
induced by uncertainties in the stellar model flux (typically 2% 
at these wavelengths). We note that we have integrated the stellar 
spectrum within each of the photometric filters before comput- 
ing the residual excess flux, i.e., the excess SED is not affected 
by a spurious color correction. 



Fig. 6. The SED of the excess emission close to the star. 
Black dots are ESO photometry, blue dots AKARI photometry. 
The red d ot refers to the near-IR interferometric measurement 
(l Absil et al. 2009), the green dot to the unresolved excess emis- 
sion at the location of the star in the Herschel PACS 70 yum im- 
age. The dashed line represents Fy oc v°-^. 

The spectral shape of the excess cannot be represented by 
a single-temperature blackbody. This could point to the pres- 
ence of hot dust that is not confined to a narrow asteroid belt, 
but is distributed over a range of distances close to the star 
(<20 AU). Alternatively, the excess may also be due to emitting 
gas. The slope of the excess emission is Fy oc yO-SiO.i^ close to 
the v"^ behavior expected for free-free emission from a spher- 
ical wind with constant outflow v elocity around a hot star, as 
derived bv iPanaeia & Fellil (Il975h . Using their equation (24), a 
reasonable stellar mass-loss rate of 10"^ M© yr"' is computed for 
Fomalhaut. Although A bsil et al.. (2009) rejected free-free emis- 
sion as the origin of the near- infrared excess based on the ab- 
sence of an emission component or variability in the Ha line, 
our determination of the slope of the SED may revive the dis- 
cussion on the origin of the excess. Further analysis is however 
beyond the scope of this paper 

4.2. Cometary debris 

The temperature of the particles in the belt, ranging from 40 K 
to 100 K in our model, together with their now known loca- 
tion, requires the presence of particles with the absorption char- 
acteristics of small grains with sizes below the blow-out size. 
Paradoxically, HST imaging showed that the dust grains in the 
belt only scatter a very smal l fraction of the stellar light they 
receive into our line of sight (iKalas et al\\2XM) . This was inter- 
preted as being due to the highly a nisotropic scattering of com- 
pact particles larger than 100 jum (iMin et a l. 2010). Very small 
compact particles, which have the required absorption properties 
to fit the temperature of the belt in our images, are incompatible 
with this very low scattering efficiency. Vice versa, the thermal 
properties of large compact 100 jum particles are inconsistent 
with the observed thermal emission in the far-infrared. 

To match both the scattered light images and our thermal in- 
frared images, we need grains that absorb and emit like small 
grains, but scatter like large particles. We suggest here that 
this combination is characteristic of fluffy aggregates, similar 
to cometary interplane tary dust particles in our Solar System 
dFraundorfet al.lll982h . Such aggregates behave in many as- 
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pects as loosely bound collections of small monomers, and thus 
can have the absorption properties of small grains, while ex- 
hibiting the s cattering anisotropy expected from large particles 
dVolten et al.li2007l) . The resulting highly anisotropic scattering 
ensures that most of the light is not scattered into our line of 
sight, giving the particles a very low apparent scattering cross 
section. 



colUsion destroys two comets with mass ^r^p^, where pc is the 
specific density of a comet. The mass-loss rate is then 



iLe Bouquin et al.l (|2009*) have resolved the stellar rotation 
axis using near-infrared interferometry. Based on the orbital di- 
rection of Fomalhaut b (i.e., counterclockwise), they argue that 
the eastern side of the belt is in fact the far side. This side is 
the bright side in scattered light, and hence the authors conclude 
that the particles are predominantly backscattering, an observa- 
tion that was elaborated upon by Min et al. (2010). However, two 
assumptions are crucial in the determination of the belt orienta- 
tion. First, Fomalhaut b must be a real planet in orbit. Second, 
the rotation axis of the planet's orbit and of the star must have 
the same direction, i.e., the planet must have a prograde orbit. If 
For nalhaut b is not a p lanet (and there is reason for skepticism, 
see iJanson et"ani2012h . the far side of the belt cannot be deter- 
mined from the available observations. If the planet is real, but 
has a retrograde orbit, the orientation is reversed, and the bright 
side would be the nearest side of the belt. In the latter case, the 
particles would be predominantly forward scattering. 

Apart from these considerations on the orientation of the 
Fomalhaut system, our "fluffy aggregate" suggestion is consis- 
tent with both forward and backward scattering, or at least not 
in conflict with either one. Large aggregates have the scattering 
properties of large grains, with a strong diffraction spike out- 
side of the observable range of scattering angles (6 e [90° - / - 
25°, 90° + i = 155°]). Most of the steUar light would be scat- 
tered forward, but in a direction close to the plane of the belt. 
Whether the belt is dominated by backward or forward scatter- 
ing at the observable scattering angles, would then depend on the 
specific scattering properties of the building blocks of the flufi'y 
aggregate. Our suggestion does exclude large compact grains, 
because these grains do not have the thermal properties needed 
to explain the far-infrared Herschel images and the SED. 

Fluffy aggregates have been invoked to explain the ob- 
served polarization p roperties of the debris disk around AU Mic 
(iGraham et al.ll200 7'). Here we used a different method to come 
to a similar conclusion on the morphology of the dust grains in 
the Fomalhaut system. 



l67T-N^v^oUPcr^^ 
3V 



dM 
dt 

W e take h-Q.\, Vcou ~ /iV orbit - 0-5 km/s (Dominik & DecinI 
l2003h . and pc = 0.6 g cm~^ (iBritt et al.ll2006h . and find that a 
population of 2.6 x lO" 10-km-sized comets, or 8.3 x 10'^ 1- 
km-sized comets are required in the ring to produce the observed 
collision rate. The Oort Cloud of the Solar System contains an 
estimated 10'^ - 1 0'^ comets which presumably came from the 
planetary system (IWeissmanlll991h . so these numbers are not 
unreasonable. The total mass of the Fomalhaut belt is about 1 10 
Earth masses, which can be compared to the approximately 30 
Earth masses required in the primordial Kuiper Belt to allow for 
the formation of Pluto and other trans-Neptunian dwarf planets 
(Stem 1996). 

Interestingly, we find the same mass in Fomalhaut's belt if 
we extrapolate the size distribution from the largest grains in 
our model (5000 pm) to the cometary regime, providing a sec- 
ond independent measurement of the same quantity. Crowding 
of planetesimals into a narrow belt through resonance capture 
by an outward-migrating planet can conceivably produce belts 
with high planetesimal content, and may have happened in our 
own Solar System's Ku iper Belt during Neptune's migration 
dHahn&Malhotralll999 ^. Large amounts of mass in the outer 
regions of disks around A-type stars like Fomalhaut might ex- 
plain the existence of planets found fa r away from the host star 
(iKalas et all 120081: iMarois et ani2008l) . as in-situ planet forma- 
tion could then have been possible. 

5. Conclusion 

The belt around Fomalhaut is a prototypical example of a steady- 
state collisional cascade. The high mass in blow-out grains evi- 
dences the extreme activity of the system. The dust particles are 
fluffy aggregates, similar in morphology to the cometary inter- 
planetary dust particles in our Solar System. This indicates that 
the planetesimals at the top of the collisional cascade are in fact 
comets. 

Acknowledgements. BV, CW, JB, WDM, KE, RH, CJ and PR acknowledge 
funding from tlie Belgian Federal Science Policy Office via the ESA-PRODEX 
office. 



4.3. Mass loss in the belt 

The amount of blow-out particles needed to reproduce the im- 
ages (Fig.O and the SED (Fig.|5]l can be directly converted into 
a mass-loss rate from the belt and, equivalently, to a dust pro- 
duction rate assuming a steady state. We find a mass-loss rate of 
2 X 10^' g yr"', corresponding to about two 10-km sized comets, 
or 2000 1-km-sized comets per day. 

We can further estimate the total number of comets in the 
ring required to generate such a high collision rate. Assuming 
that the ring ranges from ri - 133 AU to r2 - 153 AU and 
has a normalized half-height h - H/{2r), its volume is V = 
yK^I ~ '"i)- The collision rate in the belt is then given by 
A^c^coiiO"coii/(2V), where is the number of comets in the belt, 
VcoU is the relative velocity between comets, and ct^qw - 4nr^ is 
the collisional cross section of two comets with radius r^. Each 
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l20Q3h . If VcoU is the collision velocity and cr^oW is the colli- 
sional cross section, we see that the collisional time scales as 
fcoU 0^ (/(■s)'5 0"coii) We assume that VcoU is independent of 
grain size, i.e. particles with different sizes have similar orbits. 
Because ctcoW s^, we find that fcoU °^ s^^^'^"K 



PR drag pulls particles out of the ring. The inward veloc- 
ity due to PR drag at a distance r from the star is vpR oc oc 



where j3{s) is again the ratio of the radiation and 
gravitational forces on a grain of size s. At the inner radius ri 
of the ring, vpR(ri) oc /3(s). We now assume that the particle 
flux through ri into the inner disk is proportional to the mate- 
rial available in the narrow slice dr = vpR(ri) fcoii where the PR 
time scale is shorter than the collisional time scale. The particle 
flux is then proportional to F(s) oc f{s)dr - /(i) vpR(ri) fcoU 



Appendix A: Particle size and surface density 
distribution in the inner disk. 

Here we compute the size distribution of dust grains in the inner 
disk under the assumption that these grains are moving inward 
due to Poynting-Robertson (PR) drag, and that no collisions oc- 
cur in the PR stream inward of the source region. Because the 
absolute scaling (i.e. total dust mass) is determined by the fitting 
routine (see Sect. [3]l, we only consider the proportionalities of 
the different processes. 

We assume that the ring itself is in collisional equilibrium, 
so that the lifetime of particles is dominated by collisions for all 
bound particles. Even then, there will be a narrow slice close to 
the inner boundary of the ring where the time scale for a particle 
to move out of the ring by PR drag is shorter than its collisional 
destruction time scale. We assume that the particle flux seeping 
into the inner disk is proportional to the number of particles of a 
given size in that slice. Note that the width of this slice depends 
on grain size. 

In the collisional ring itself, we assume that the size distribu- 
tion is given by f(s) oc s", where a = -3.5 in case of collisional 
equilibrium. First, we estimate the collisional time scale for par- 
ticles of size s. The most important collision partners are grains 
with sizes near s, in a size range ds oc s dPominik & DecinI 



Mass conservation in the flow through the inner disk then 
gives the surface density E(i) - where we now have 

to consider that vrr is no longer evaluated at the inner ra- 
dius of the ring, but changes proportional to ^ . We find that 
l,(s) oc = s^^ ■ For an undisturbed PR flow, the surface den- 

sity depends on the grain size, but not on the distance to the star 
Based on this reasoning, we use a size distribution f{s) oc s^^ in 
the inner disk. Note, however, that in reality, the size distribution 
may be modified by interactions with a planetary system. Also 
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other sources of dust, such as evaporating comets or an unre- 
solved inner belt, would give rise to different size distributions. 
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